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Abstract. We report on the analysis of ~ 22,000 M dwarfs using a statistical parallax 
method. This technique employs a maximum-likelihood formulation to simultaneously 
solve for the absolute magnitude, velocity ellipsoid parameters and reflex solar motion 
of a homogeneous stellar sample, and has previously been applied to Galactic RR Lyrae 
and Cepheid populations and to the Palomar/Michigan State University (PMSU) survey 
of nearby low-mass stars. We analyze subsamples of the most recent spectroscopic cat- 
alog of M dwarfs in the Sloan Digital Sky Survey (SDSS) to determine absolute magni- 
tudes and kinematic properties as a function of spectral type, color, chromospheric ac- 
tivity and metallicity. We find new, independent spectral type-absolute magnitude rela- 
tions, and color-absolute magnitude relations in the SDSS filters, and compare to those 
found from other methods. Active stars have brighter absolute magnitudes and lower 
metallicity stars have fainter absolute magnitudes for stars of type M0-M4. Our kine- 
matic analysis confirms previous results for the solar motion and velocity dispersions, 
with more distant stars possessing larger peculiar motions, and chromospherically ac- 
tive (younger) stars having smaller velocity dispersions than their inactive counterparts. 
We find some evidence for systematic diff'erences in the mean U and W velocities of 
samples subdivided by color 



1. Introduction 



M dw arfs are the most numerous stellar population in the Milky Wav (IBochanski et al. 
2010I) . Surveys such as the Sloan Digita l Sky Survey (SDSS ; lYork et alJ l2000l) and 



the Two Micron All Sky Survey (2MASS: ISkrutskie et aDl2006h have led to photomet- 
ric and (in the case of SDSS) spectroscopic catalogs contai ning larg e numbers of low 
mass stars. The largest spectroscopic database of M dwarfs dWest et al. 2010) compiles 
spectral types, colors, proper motions, radial velocities and chromospheric activity es- 
timates (as traced by Balmer series emission) for more than 70,000 stars from the most 
recent SDSS data release. 

In order to make widespread use of these new measurements of the low-mass 
stellar population, with a focus on Galactic structure and kinematics, it is necessary 
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to have good distance estimates. Since only very nearby M dwarfs have measured 
trigonometric parallaxes, photometric and spectroscopic parallax relations have typi- 
cally been employed to obtain distances based on a star's color or spectral type re- 



spect ively (iHawlev et al.ll2002l : IWest et al.ll2005l : lKraus & Hillenbrandll2007l : IJuric et al. 
l2008h . Because SDSS photometry saturates at m ~ 15, it is particularly difficult to an- 
chor the photometric parallax relations in the SDSS filters with measured trigonometric 
parallax stars. 

The classical statistical parallax method is a way to determine the absolute mag- 
nitude of a homogeneous set of stars. Statistical parallax analysis seeks to determine 
the distance scale that provides the best match between the measured proper motions 
and radial velocities of a given stellar population, returning an estimate of the average 
absolute magnitude of the population, and the kinematic properties including the reflex 
solar motion and the velocity ellipsoid (velocity dispersions along three principal axes). 
Previo u s discus s ions o f the statistical parallax method can be found in IClube & Jones 
(Il97lh . iMurravl (Il983h and |Popowski & Gould (fT 998i). Our particular formulation has 
been used to study RR Ly raes (Hawlev et al. 198 61: IStrugnell et all 1986 : Lavden et al. 



19961 : iFemlev et al.lll998l : IPopowski & Gouldl 119981 ). Cepheids (IWilson et al.lll99lh . 



and the nearby low-mass stars from the PMSU survey (IHawlev et al.ll 19961) . 



2. Observations 

Accurate photometry, proper motions and radial velocities are required for statistical 
parallax analysi s. Th e data were obtained from the latest SDSS data release (DRV; 
IXbazajian et all |2009|) which contains photom etry over nearly 10,000 sq. deg. down 
to r ~22 in five filters {ugriz, iFukugita et al.lll996l) . When sky conditions at Apache 
Point Observatory were not photometric, the SDSS operated in a spectroscopic mode. 
Two fiber-fed spectrographs collected 640 spectra simultaneously with a typical total 
exposure time of ~ 45 minutes. These medium-resolution {R ~ 2, 000) spectra cover 
the entire optical bandpass (3800 - 9200 A). Approximately 460,000 stellar spectra are 
available in the DRV data release, obtained both as targeted and serendipitous observa- 
tions. 



^The absolute astrometric precision of SDSS is < 0. 1" in each coordinate (iPier et al 

2003h. Proper rnotions for SDSS objects were found by matching to the USNO-B sur- 
vey (iMunn et al.ll2004 ). The proper motions have a baseline of ~ 50 years, and a typical 
precision of ~ 3 mas yr~' in right ascension and declination. For a typical thin disk star 
with a transverse velocity of ~ 10 km s this limit corresponds to a limiting distance 
of ~ VOO pc. We therefore restricted the sample to stars closer than VOOpc as calcu- 
lated from the color-magnitude relation of iBochanski et al.l (l2010l) . to avoid biasing the 
solution by including only high velocity stars at larger distances. 



The lWest etall (l2010l) SDSS M dwarf sample was obtained by first selecting all 



stars in the DRV spectroscopic cat alog with colors r - / > 0.42 and i - z > 0.24, before 
correcting for Galactic reddening (jSchlegel et al."l998'). Each spectrum was processed 
through the Hammer IDL package (Covey et al. 200V) which provides a spectral type 
estimate and measures a number of spectr al features, including Ha equivalent width and 
various TiO and CaH bandhead strengths. Each spectrum was also visually inspected 
and non-M dwarf contaminants were culled from the catalog. The radial velocities 
were measured by cros s-correlating each spectrum with the appropriate template from 



Bochanski et all (I200VI) . with a precision of ~ V km s . To obtain the sample used in 
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our statistical parallax analysis, w e required precise photo metry and accurate proper 
motions and radial velocities. See Bochans ki et al. (I2010h for information on photo- 
metric flag cuts, and Dhital et al. (2010) for our kinematic quality flags. After these 
quality cuts, our sample contained 40,963 stars. The 700pc distance cut further reduced 
the final statistical parallax sample to 22,542 stars. 



3. Method and Sample Subdivision 

We employed the m aximum-hkelih ood formulation of classical statistical parallax anal- 
ysis as presented by iMurrayl ( 19831) . Briefly, the velocity distribution of a homogeneous 



stellar population is modeled with nine kinematic parameters, including the three or- 
thogonal velocities of the reflex solar motion, and the three directions and three dis- 
persions of the velocity ellipsoid, which describes the random and peculiar velocities 
of the population. Two distance scale parameters allow the absolute magnitude and 
its dispersion to vary; these are used to transform the observed proper motions into 
transverse velocities. There are thus eleven parameters used in the model. In order to 
achieve well-determined solutions, we fixed the dispersion in the absolute magn itude 
to be (Tm - 0.4 which is the observed value for low-mass stars i Bochanski et al.l l2010). 

The observational data needed are the position, apparent magnitude (corrected for 
Galactic extinction), proper motions and radial velocity for each stai" in the sample. We 
solve for the model parameter s simulta neously by maximizing the likelihood using geo- 
metric simplex optimization dNelder & Meadlll965l : lDanielslll978h . Uncertainties in the 
parameters are determined by numerical computation of the derivative of each param- 
eter individually, while keeping all other parameters fixe d. The maxirnum-lik elihood 
equations and simplex technique are described in detail in lHawley et al.l ( 19861) . 



Our spectroscopic sample is much larger than was available for previous statisti- 
cal parallax studies, which typically contained observations of a few hundred or less 
objects. Thus, we were able to divide the sample into smaller samples subdivided by 
color, spectral type, magnetic activity (as traced by Ha), and metallicity (using the ^ 
parameter of Lepi ne et al . 2007). Each subsample of M dwarfs was analyzed using the 
following prescription. Ten runs were calculated for each dataset. The initial abso- 
lute magnitude e stimates were derived from the M^, r - z color-magnitude relation of 



Bochanski et all (12010) . The absolute magnitude estimate was updated after each run. 



with the output of the previous run being used as input for the next. For a given run, 
the simplex optimization iterated 5,000 times. Thus, for each subsample, 50,000 iter- 
ations were computed. This ensured that the simplex operator had sufficient freedom 
to explore parameter space and converged to the best (maximum likelihood) answer. 
Typically, convergence was obtained after 25,000 iterations. 



4. Results: Absolute Magnitudes 

Traditionally, colors or spectral types have been employed to estimate absolute magni- 
tude when trigonometric parallaxes were not available. These are referred to as pho- 
tometric and spectroscopic paraUaxes, respectively. In Figure 1 (left panel) we show 
the spectroscopic parallax relation for M dwarfs and cornpare it to previous studies 
(iHawlev et al.ll2002l : IWest et al.ll2005l : lKraus & HiUenbrandll2007l) . The open circles are 
nearby low-mass stars with accurate trigonometric parallaxes, tabulated in iBochanski 
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( 2008h . Our results highlight the large dispersion in M,. as a function of spectral type. 
The typical spread is ~1 mag, increasing to ~3 mag near type M4, much larger than the 
spread for a given r-z color bin (see Figure 1, right panel). This suggests that spectral 
type is a poor tracer of r-band absolute magnitude for M dwarfs, since a single spectral 
type can encompass such a range of values. 

Figure 1 (left panel) also displays disagreement between our average statistical 
parallax result (solid black line, filled circles) and the nearby star sample (open circles) 
at types later than about M4. This is due to a systematic color difference within a 
given spectral type bin between the two samples. That is, the nearby star sample is 
systematically redder than the SDSS sample at a given spectral type, possibly due to 
the small number of nearby stars at late types. Thus, the SDSS stars have an absolute 
magnitude (at a given spectral type) that is appropriate for the bluer, more luminous 
stars of that type. 

Figure 1 (right panel) give s the statistical parallax Mr, r-z relation, along with 
the re sults of previous studies (Hawley et al.ll2002l : IWest et al.ll2008l : iBochanski et al. 



120101 '). Overplotted with open circles are the absolute magnitudes and colors from the 
nearby star sample. Note that the dispersion in absolute magnitude as a function of 
color is significantly smaller than that shown in Figure 1 (left panel) indicating that 
the r-z color is a much better tracer of r-band absolute magnitude. Furthermore, the 
discrepancy in mean color for a given spectral type between the nearby star sample and 
the SDSS stars is not evident. We reiterate that colors, rather than spectral types, are 
preferred for absolute magnitude and distance estimation for M dwarfs. 
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Figure 1. Left Panel - Mr vs. spectral type. The symbols are described in the text. 
Right Panel - M,- vs. r-z, same symbols as in left panel. 

Due to the large spread in Mr as a function of spectral type, we concentrate on 
the color results shown in Figure 1 (right panel) to investigate magnetic activity effects. 
Our results are shown with the black, red and blue lines. There are a few notable trends. 
First, the bifurcated main sequence at blue colors (earlier spectral types) indicates that 
active, early-type M dwarfs are about ~ 1 mag brighter in Mr than inactive stars of the 
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same color. A similar effect was observed in the P MSU sample, where active M dwarfs 
were brighter in V than inactive stars (Fig ure 4 of iHawlev et alJll996l) . This difference 
could be explained if the active stars have larger radii, as suggested b y recent observa- 
tions of eclipsing binarie s ( Lopez-Mo rales 200?!: [Morales et aU2010l) and active, single 
stars (^Berger et al.''2006'; Morales et al. 2008). Direct interferometric measurements of 
stellar radii will be necessary to independently quantify the amount the radius changes 
for an individual star, since our statistical analysis is performed on subsamples of stars 
that span a small range of temperature. 

At blue colors {r - z < 2), the inactive and total sample absolute magnitudes fall 
below the mean locus of nearby trigonometric parallax stars. Since the nearby stars are 
composed of a mix of active and inactive stars, activity is likely not the only important 
effect on the luminosity. 

We next investigated the effects of metallicity on Mr and attempted to isolate them 
from those traced by chromospheric activity. Metallicity can alter a star's effective 
temperature and luminosity, manifesting as a change in absolute magnitude. At a given 
color (or spectral type), stars with lower metallicity typically exhibit fainter absolute 
m agnitudes. We exarn ined the effects of metallicity using the ^ parameter, as defined 
bv lLepine et"aD(l2007h . ^ is a metallicity proxy that can be employed as a rough tracer 
of [Fe/H], wi th 4" = 1 correspon ding to solar metallicity and ^ = 0.4 corresponding to 
[Fe/H] ~ -1 (IWoolf et al.ll2009h . We note that the [Fe/H], ^ relation is only calibrated 
for early spectral types (~ M0-M3) and has significant uncertainty near solar metallicity. 

In Figure 2, we plot the absolute magnitudes for two values of ^ as a function of 
spectral type (left panel) and r - z color (right panel) for active and inactive stars. The 
uncertainty in each bin is cr « 0.4 mags which was fixed in the solution (see above). 
Figure 2 demonstrates that for stars with similar chromospheric properties, the lower 
metallicity ones have fainter absolute magnitudes at the same color or spectral type, as 
expected. Active stars at the same metallicity are brighter than their inactive counter- 
parts, but both metallicity and activity are important for determining the luminosity of 
an individual star. At a fixed ^, the difference in Mr is consistent with Figure 1, with 
early-type active stars being ~ 1 mag brighter, and the difference diminishing at later 
spectral types. The total SDSS sample is overplotted in each panel with a solid black 
line. The early-type stars, which are seen at larger distances due to SDSS magnitude 
limits, fall near the lower metallicity, inactive loci. The later type stars, located closer 
to the Sun, are consistent with solar metallicities and active stars. This suggests that the 
low-mass dwarfs are tracing a met alhcity gradient si milar to the one observed in SDSS 
observations of higher-mass stars (llvezic et al.ll2008h . and also explains why the SDSS 
sample falls below the locus of nearby stars at early types in Figure 1. 

We also note that the activity-metallicity loci appear to converge near type M5 
(r - z ~ 2.8). This behavior is not well explained, but may be linked to the transition 
between a partially and fully convective stellar interior that occurs near that spectral 
type/color. Perhaps this transition, which alters the efficiency of energy transport in the 
star, also regulates the luminosity at the surface. 



5. Results: Kinematics 



We measure kinematics in a WW coordinate system with U increasing towards the 
Galactic center, V increasing in the directio n of solar motion, and W increasing verti- 
cally upward (as in lDehnen & Binneylll998h . 
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Figure 2. Left Panel - M,- vs. spectral type for two values of ^, a metallicity proxy. 
The total sample is overplotted in the solid black line. Right Panel - The Mr vs. r-z 
results show the same features as in the spectral type panel. 



The statistical parallax method gives the reflex solar motion with respect to the 
mean velocity of the sample. If samples possess different mean velocities, this will be 
manifested as a change in the solar reflex motions that are returned. In the U and W 
directions, the solar reflex motion reflects the Sun's peculiar motion, while in the V 
direction it is a combination of the Sun's peculiar motion and the asymmetric drift at 
the solar circle. This increases the mean V velocity of a sample of typical disk-age M 
dwarfs, giving a V reflex motion that is larger than one measured from a population of 
young stars. Our results are shown in Figure 3 (left panel). While the U and W veloc- 
ities remain relatively constant with spectral type, the V reflex motion is significantly 
larger at bluer colors. These stars are observed at larger distances, and also have larger 
velocity dispersions (as seen in Figure 3, right panel). For r - z > 2, the V velocity 
of the Sun remains relatively constant at 20 km s~^ This value compares favorably 
to the V velocities previously measured for M dwarfs dHawley et alJll996l : iFuchs et all 
12009) . At later types, where the separation in age is most extreme between the active 
and inactive populations, the active stars in the sample show decreasing V velocities, 
while the inactive stars have increasing V velocities as expected. 

While the U and W velocities should be constant with color (or spectral type), we 
find that there is significant structure in both distributions. The W velocity distribution 
in Figure 3 (left panel) appears to peak near r - z ~ 2.3 (type M4). This structure indi- 
cates the mean vertical motions of the M dwarf subsamples are varying, with both bluer 
and redder M dwarfs having smaller mean W velocities. The W velocities of active and 
inactive stars are not significantly different. Meanwhile, the U velocity distribution 
shows a different behavior, exhibiting a rise toward later types which begins at a bluer 
color in the active stars. 

The active M dwarfs are k nown to forrn a kinematically colder population than 
their inactive counterparts (e.g.. IWielenlll977l) . as evidenced by smaller velocity dis- 
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Figure 3. Left Panel - Solar peculiar motion as a function of r - z color. There 
is unexpected structure in the U and W distributions. Right Panel - UVW velocity 
dispersions vs. color for active (red triangles), inactive (blue squares) and all stars 
(filled black circles). Active stars possess smaller dispersions at red colors. 



persions. This is usually interpreted as an age effect, where younger M dwarfs have 
experienced less dynamical heating. Activity in low-mass stars has been shown to de- 
pend on both spectral type and age; i nactive, late-typ e M dwarfs are usually much older 
than inactive early-type M dwarfs ( West et aP 2008). Figure 3 (right panel) shows our 
velocity dispersion results. As expected, the dispersions of late-type (younger) active 
stars are smaller than for the late-type (older) inactive stars. The overall decline in ve- 
locity dispersion toward redder colors reflects a similar effect, since the early type stars 
are sampled further away, and are therefore older and have higher dispersions than the 
closer, later-type stars. 
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